The effect of disorder on transverse domain wall dynamics in magnetic nanostrips 
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We study the effect of disorder on the dynamics of a transverse domain wall in ferromagnetic 
nanostrips, driven either by magnetic fields or spin-polarized currents, by performing a large en- 
semble of GPU-accelerated micromagnetic simulations. Disorder is modeled by including small, 
randomly distributed non-magnetic voids in the system. Studying the domain wall velocity as a 
function of the applied field and current density reveals fundamental differences in the domain wall 
dynamics induced by these two modes of driving: For the field-driven case, we identify two different 
domain wall pinning mechanisms, operating below and above the Walker breakdown, respectively, 
whereas for the current-driven case pinning is absent above the Walker breakdown. Increasing the 
disorder strength induces a larger Walker breakdown field and current, and leads to decreased and 
increased domain wall velocities at the breakdown field and current, respectively. Furthermore, for 
adiabatic spin transfer torque, the intrinsic pinning mechanism is found to be suppressed by disor- 
der. We explain these findings within the one-dimensional model in terms of an effective damping 
parameter a* increasing with the disorder strength. 

PACS numbers: 75.78.Fg, 72.25.Ba, 75.78.Cd 



Domain wall (DW) dynamics in nanoscale ferromag- 
netic wires and strips driven by magnetic fields or spin- 
polarized currents is a subject of major technological im- 
portance for the operation of potential future nanoscale 
magnetic memory [l|,0] and logic 0] devices. In these de- 
vices information is typically stored as magnetic domains 
along a nanostrip/wire and is processed by DW motion. 
For the reliable operation of such devices it is of funda- 
mental importance to understand and control the effect 
of imperfections or disorder on the DW dynamics, nec- 
essarily present in any realistic samples, e.g. in the form 
of thickness fluctuations and grain structure of the sam- 
ple, or various impurities and defects in the material. At 
the same time, such systems constitute a low-dimensional 
limit of the general problem of driven elastic manifolds 
in a random potential Q. 

While the crucial importance of disorder for the dy- 
namics of higher-dimensional DWs is well established, 
resulting in phenomena such as the Barkhausen effect 
0, majority of studies of DW motion in systems with 
nanostrip/wire geometry neglect disorder effects. This 
applies to both theoretical studies and interpretations of 
experimental results. Some exceptions include studies 
demonstrating enhanced DW propagation due to rough- 
ness of the edges of the strip [y, |7| ■ Recently also the ef- 
fect of spatially varying saturation magnetization M s on 
the dynamics of vortex walls was studied, resulting in an 
effective damping increasing with the disorder strength 
0. Similar spatially distributed disorder has also been 
studied in a simplified, line-based model of a transverse 
DW @, H(ij |. Experimental studies of DW dynamics in 
wires have revealed its stochastic nature in the case of 



presence of disorder in the samples, in combination with 
thermal effects. For longer current pulses, the resulting 
average DW velocities have been shown to be quite low 
[IH, likely due to pinning effects induced by structural 
disorder. Dynamical pinning effects have also been ob- 
served in experiments of field-driven vortex wall dynam- 
ics HI 14 1 . However, despite of these advances, many 



short current pulses [111 ], and has been attributed to the 



details of the disorder effects on DW dynamics in nanos- 
tructures remain to be clarified. 

In this Letter, we consider by micromagnetic simu- 
lations the effect of disorder on the field and current- 
driven dynamics of a transverse DW in a narrow and 
thin Permalloy strip. Disorder is modelled by including 
randomly positioned small non-magnetic regions (voids) 
in the system. Our results show that the field and 
current-driven DW dynamics exhibit remarkable differ- 
ences which are only revealed in the presence of disor- 
der. In particular, we identify two fundamentally dif- 
ferent DW pinning mechanisms acting in a field-driven 
system, operating below and above the Walker break- 
down field, respectively, with the latter mechanism being 
absent in the current-driven case. Also the Walker break- 
down itself is affected by the presence of disorder, such 
that it is shifted to larger field and current values with 
increasing disorder strength. At the same time the DW 
velocities at the breakdown field and current get smaller 
and larger, respectively. Furthermore, for adiabatic spin 
transfer torque, the intrinsic pinning mechanism is found 
to be suppressed by disorder. These findings emphasize 
the importance of understanding the interplay between 
disorder, the DW structure and the properties of the ex- 
ternal driving force, and are shown to be related to an 
effective damping parameter a* increasing with the dis- 
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FIG. 1. (Color online) The average velocity v m of the moving 
DWs (main figure) and v emp = (1 — Ppm)v m (inset) as a func- 
tion of H sx t and a. Error bars correspond to the standard 
deviation of v m - The pinning probabilities P p i n during the 
20 ns simulation (bottom panel) exhibit large values for large 
Hext due to the core pinning mechanism. 



order strength. 

We perform a large ensemble of micromagnetic sim- 
ulations with the GPU-based micromagnetic simulator 
MuMax 15[ , making it possible to obtain large statistics 
for averaging over the disorder realizations. To study the 
time evolution of the magnetization M(r, t) with an am- 
plitude M s , we solve the Landau-Lifshitz (LL) equation 
with the spin-transfer torque terms [lij] . 
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where H e ff is the effective magnetic field (with contri- 
butions from the external, exchange and demagnetiza- 
tion fields), 7 is the gyromagnetic ratio, a is the Gilbert 
damping constant, £ is the degree of non-adiabaticity, 
j is the current density, and bj — P//s/(eM s (l + £ 2 )), 
with P the polarization, (is the Bohr magneton and e 
the electron charge. 

We consider Permalloy strips of width w — 100 urn and 
thickness 10 nm, such that the stable DW structure is a 
head-to-head V-shaped symmetric transverse wall, sepa- 
rating in-plane domains pointing along the strip axis [17| . 
The used material parameters are those of Permalloy, i.e. 
M s = 860 x 10 3 A/m and a — 0.02, and no anisotropy 
fields are included in Eq. (p}. To clearly see the effect 
of quenched disorder on the DW dynamics, we set the 
temperature T = 0. The system is discretized by con- 
sidering N cells of size 3.125 x 3.125 x 10 nm 3 . Upon 



application of an external magnetic field H ex t = H ext x. 
along the strip axis in the absence of disorder, the DW is 
displaced along the strip. If the field is below the Walker 
breakdown field Hw, the DW essentially keeps its equi- 
librium structure during the propagation, with a small 
out-of-plane component close to the tip of the V-shape, 
and a velocity roughly linearly proportional to the ap- 
plied field. Above Hw a- n antivortex is nucleated at the 
tip of the V-shape. It then propagates across the strip 
width, reversing the polarity of the DW magnetization. 
This process is repeated such that the DW polarity oscil- 
lates back and forth, dramatically decreasing the average 
DW velocity 0. 

With disorder included in the form of randomly posi- 
tioned non-magnetic voids of linear size 3.125 nm with 
varying densities a within a strip of length L — 3.2 /mi, 
the DW can get pinned even for non-zero applied fields 
[2fj| . This makes measurement and even definition of the 
DW velocity a non-trivial task. Thus, in what follows we 
consider both the "conditional velocities" v m of the mov- 
ing DWs, conditioned on the fact that the DWs will not 
get pinned during the time interval At = 20 ns we con- 
sider in the simulations (i.e. the DW will either reach the 
end of the strip or it is still moving after At = 20 ns) [lj| , 
and the probability P V i n for the DW to get pinned during 
At. These are computed by averaging over 50 disorder 
realizations for each H ext and a. Notice that here we con- 
sider a T = system, such that a pinned DW cannot de- 
pin. An alternative measure of the DW velocity (which is 
likely to be closer to typical experimental measurements 



where T > 0) is given by v e 



{l-P P in)v m - In general, 



Ppin will increase with the observation (time and length) 
scale, thus making also v exp a scale-dependent quantity. 

Fig. Q] shows the resulting average velocities v m of the 
moving DWs as a function of H ext and a. The presence 
of voids induces a finite depinning field Hd ep {&) increas- 
ing with a. For H ext > Hde P (&), v m first increases until 
a maximum velocity is reached at H ext — H\v(&), and 
then starts to decrease again. The position Hw{v) of 
this maximum, corresponding to the Walker breakdown, 
is shifted towards larger field values as a is increased, 
and the corresponding maximum velocity v m (Hw(o')) 
decreases with a. The error bars in Fig. [T] correspond to 
the standard deviation of v m , and indicate that the dy- 
namics of moving DWs has a stochastic nature due to the 
random disorder. Notice in particular that the pinning 
probability P p i n exhibits a non-monotonic dependence on 
Hext, with strong pinning for both small and large H extl 
while for intermediate applied fields (corresponding to 
large values of v m ) pinning is less likely. The maximum 
value of v exp (inset of Fig. [T]) exhibits a strong depen- 
dence on cr, and depends also on the observation scale 
via Ppi n (not shown). For large H ext , P p i n is close to 1 
for A< = 20 ns, and consequently v exp is essentially zero. 
Similar pinning effects for large applied fields have been 
observed experimentally for vortex walls [HI , . 
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FIG. 2. (Color online) Examples of the spatial distribution of 
the contributions of the applied field H ex t = 5 mT (top) and 
current density j ext = 20 x 10 12 A/m 2 with £ = (middle) 
to dlSA/dt in Eq. {TJ, corresponding to the magnetization 
configuration shown in the bottom panel, exhibiting an an- 
tivortex in the middle of the strip. dM./dt is given in units 
of M„/s. The randomly positioned voids with a = 3125^im~ 2 
are shown as grey dots. 



To gain insight on the mechanisms behind this behav- 
ior, we consider snapshots of the DW configurations and 
the various contributions to <9M / dt in Eq. (QJ . For small 
Hext, we find that the overall DW structure is preserved, 
with the disorder inducing only minor distortions. If 
the DW gets pinned, this happens by a collective ac- 
tion of several voids. This mechanism is known as col- 
lective pinning, and it is responsible for the non-zero 
depinning field Hd ep < Hw(cr). Remarkably, we identify 
a fundamentally different pinning mechanism for large 
fields, H ext > Hw(&): In this regime, an antivortex is 
able to propagate to the interior of the strip, resulting 
in pinned DW configurations (occurring with probabil- 
ity Ppin) with the antivortex core positioned exactly on 
top of a void or a local void structure. We refer to this 
mechanism as core pinning, and attribute it to the fact 
that the energy of the system can be significantly lower 
when the antivortex core or part of it - involving large 
magnetization gradients and out-of-plane magnetization 
- is placed in a non-magnetic region (or more generally, 
in a region with low M s ). In the field-driven case the DW 
is susceptible to get pinned by this mechanism because 
the Zeeman torque is relatively small in magnitude and 



does not directly displace the DW (top panel of Fig. [2]); 
Instead, the small out-of-plane magnetization due to the 
Zeeman torque induces demagnetizing fields, which act 
to move the DW. Such an indirect driving mechanism is 
sensitive to the perturbations due to disorder, leading to 
several effects, including a-dcpcndent H c i ep and Hw, and 
in particular the core pinning mechanism for high H ext . 

We proceed to contrast these results with the current- 
driven case, by applying a current density j = —j ex t* 
with P — 0.5 along strips of length L=6A /im. We first 
consider perfect adiabaticity (£ = 0, top panel of Fig. 
131). Due to intrinsic pinning 2l|, there is a non-zero de- 
pinning current jdep,int in the absence of disorder, above 
which DW motion involves repeated polarity transfor- 
mations mediated by antivortex propagation across the 
strip width. Adding disorder with the same procedure 
as above reveals two intriguing observations: First, it 
appears that the DW is able to move even for currents 
slightly below jdep,int ■ This surprising finding can be ex- 
plained by noticing that the intrinsic pinning mechanism 
is due to the ability of the DW to deform in such a way 
that the torques due to interactions within the DW (i.e. 
the effective field) exactly counterbalance the adiabatic 
spin-transfer torque 21] . However, the presence of dis- 
order induces additional DW deformations and imposes 
constraints on the ability of the DW to counteract the 
current-induced torques, leading to non-zero values for 
both v m and 1 - P pin for j ext somewhat below jd ep ,int- 
Notice that while v exp (inset of the top panel in Fig. |3j) 
exhibits non-linear field dependence reminiscent of typi- 
cal creep motion for small fields, we are considering here 
aT = system in which a pinned DW cannot depin due 
to the absence of thermal fluctuations [22[ . 

The second observation is that for larger j ext , core pin- 
ning is absent. Even if for j ext > jw{o~) the antivor- 
tex core is constantly moving back and forth across the 
strip width, it never gets pinned by the voids, strongly 
contrasting with the field-driven case. To explain this 
observation, we consider the spatial distribution of the 
current-induced contribution to dM./dt (middle panel of 
Fig. [5]), and find that the current acts directly (in con- 
trast to the indirect mechanism in the field-driven case) 
and strongly on the antivortex core where the magnetiza- 
tion gradients are large, facilitating its propagation along 
the strip across the energy barriers due to the voids. This 
is also directly visible in the the LL equation (Eq. ((TJ)), 
where the current acts on the gradient of M rather than 
on M itself. 

Finally we consider the role of the non-adiabatic spin- 
transfer torque (bottom panel of Fig. [31 where the 
£ = 0.04 case is shown) on the DW dynamics. For £ > 
and a = 0, there is no intrinsic pinning, and the DW 
propagates preserving its internal structure with a finite 
velocity linearly proportional to the current density j ext 
up to a Walker breakdown current jw ■ For j ext > jw , an 
antivortex is again nucleated and propagates across the 
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FIG. 3. (Color online) The average velocity v m of the moving 
DWs as a function of jext and a, for £ = (top) and £ = 0.04 
(bottom). Error bars correspond to the standard deviation 
of v m . The pinning probabilities P P i„ during the 20 ns simu- 
lation highlight the absence of core pinning for large current 
densities. The insets show v e x P = (1 — P p in)v m for £ = (top 
panel), and the effective a* (a) for various £ (bottom panel), 
respectively. 



strip width reversing the polarity of the DW magnetiza- 
tion, and decreasing the average DW velocity. For larger 
jext, the velocity again increases with j ext . Adding dis- 
order induces a finite depinning threshold jd ep (cr), and 
pushes the local maximum of v m or the Walker break- 
down to higher j ext . At the same time, v m at jw(&) 
increases with a. Thus, the voids are able to inhibit the 
antivortex entering the strip, enhancing the DW prop- 
agation and structural stability for intermediate current 
densities, jwip = 0) < j ext < jwip > 0). This effect 
arises as the antivortex core is pushed across the strip 
width by the effective field terms in Eq. (TTJ) (notice that 
the effect of the current is symmetric such that no an- 
tivortex displacement along the y direction arises directly 
due to the current, see the middle panel of Fig. [5]), a 
mechanism sensitive to the disturbances due to disorder. 
Again, there is no core pinning for j ext > jw{&), for the 
same reason as in the adiabatic (£ = 0) case. 



For jdep(o-) < jext < jw(cr), v m depends linearly 
on jext, and by extrapolating linear fits to the data to 
jext = all the lines cross at v m — (not shown). Thus, 
we estimate effective values of the damping parameter 
from the slopes of these linear fits [8j, as within one- 
dimensional models [23] v m oc {j3/a)j e xt for jext < jw, 



with j3 = £/(l +£ 2 ). Our simulations (inset of the lower 
panel of Fig. [3]) with different £ indicate that the data 
can be interpreted in terms of an effective a* increas- 
ing with cr [§[. Also an effective M* = (1 - aLw/N)M s 
emerges naturally. Thus we can explain our results with 
the one-dimensional model in terms of cr-dcpcndent effec- 
tive parameters: For instance, jw(&) — iTTj(Mg &\N y — 
N x \)*a* /(gn B P\P - a*\), with A the DW width and 
N x and N y the demagnetizing factors, and jdep,int{o') = 
jw(&,£, = 0) 23]. Using the expression for jw and 
the values of a* to estimate C* = (AM^\N y - N x \)*, 
the scaling of jdep.int with a can be reproduced re- 
markably well, see Table [IJ A similar analysis in the 
field-driven case, with H w = 2ira*(M s \N y — N x \)* and 
Vm(Hw) = (7 A* /a*)Hw [13 , reproduces the observed 
scaling of both Hw and v m (Hw) with a (Table IJ). No- 
tice that in our case v m (Hw) depends on a through 
the CT-dependent effective parameters, while for systems 
with only edge roughness v m (Hw) is independent of the 
amount of edge roughness [ff. 

To summarize, we have presented a detailed analysis 
of the effect of disorder on the field and current-driven 
transverse DW dynamics in a narrow and thin Permalloy 
nanostrip. We have identified two fundamentally differ- 
ent pinning mechanisms, acting in different regimes of the 
DW propagation. The observation that there is no core 
pinning in the current-driven case whereas it dominates 
the field driven dynamics for large fields highlights the 
different nature of the field and current drive in a way 
that can be observed only in the presence of disorder. In 
general, we have seen that the pinning mechanisms oper- 
ating will depend on the details of the DW structure, and 
thus we expect that the core pinning mechanism is absent 
for systems with high perpendicular magnetocrystalline 
anisotropy as there is no (anti)vortex core that could get 
pinned, but it could play a role in the dynamics of vortex 
walls occurring in wider soft strips Q, possibly also for 
small applied fields. If only edge roughness is present, no 
core pinning should occur. Experiments should be per- 
formed to systematically study the scale dependence of 
P p in and v exp . Finally, we point out that the observation 
that disorder tends to stabilize the DW internal structure 
and increase the maximum DW velocity by suppressing 
the Walker breakdown in the current-driven case suggests 
that it could be desirable to deliberately engineer disor- 
der in the system, for instance to replace notches to pin 
the DW in various technological applications [24j ■ 

Acknowledgments. Stefano Zapperi is thanked for 
numerous interesting discussions on DW dynamics and 
disorder, and Mikko J. Alava for useful comments on the 



5 



a [/mi 2 ] 


a* 


C* [A 2 /m] 


3dep,int [ A / m 1 


■sim r A / 2i 

Jde P ,int [A/m J 


H% ed [mT] 


H^ m [mT] 


v^ ed {Hw) [m/s] 


C(%) [m/s] 





0.0200 


2.92 xlO" 10 


14 xlO 12 


15 xlO 12 


2.75 


2.75 


457 


457 


1562.5 


0.0221 


2.52 xlO" 10 


12.1 xlO 12 


13 xlO 12 


3.05 


3.0 


398 


437 


3125 


0.0238 


2.45 xlO -10 


11.7 xlO 12 


12 xlO 12 


3.25 


3.25 


389 


419 


4687.5 


0.0258 


2,36 xlO" 10 


11.3 xlO 12 


11.5 xlO 12 


3.52 


3.25 


377 


403 


6250 


0.0283 


2.28 xlO" 10 


10.9 xlO 12 


11 xlO 12 


3.78 


3.5 


368 


387 



TABLE I. Predictions for idep.int, Hw and v m (Hw) from the one-dimensional model in terms of er-dependent effective a* and 
C* = (AMs\N y — N x \)* , compared with the simulated values. C* is estimated by fitting the expression for jw (see text) to 
the data in the bottom panel of Fig. [3] 
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